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ABSTRACT 

Hydraulic ram concerns catastrophic failure of liquid 
filled tanks impacted by high speed projectiles. It is 
divided into four phases: entry/shock, pressure field, 

cavity, and free surface. The entry wall phase was investi- 
gated in detail. 

A formula to predict energy losses of a bullet penetrating 
an aluminum plate is presented, and experimental data were 
obatined and compared. 

A fiber optic transducer was designed to measure generated 
pressures in the shock and cavity phases. 
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I . INTRODUCTION 



Since the days aircraft first began to fly in a hostile 
environment, the achievement of non-nuclear aircraft sur- 
vivability has been the subject of extensive research. The 
ability of today’s expensive air-to-air and air-to-ground 
weapons delivery systems to withstand hits by small arms 
fire and by warhead fragments is of paramount importance in 
this era of decreasing defense budgets. 

An aircraft component which is easily hit and most 
susceptible to damage is the fuel tank. A damaged fuel cell 
can lead to aircraft loss in several ways, some of which are 
fire, explosion, fuel starvation, or fuel cell failure. 

Catastrophic failure of a fuel tank due to impact by a 
high speed projectile is known as hydraulic ram. This is 
actually a misnomer, since the failure mechanism consists of 
four distinct phases rather than one quasi-static compression 
type of state. [Ref. 1] 

The first phase, the entry/shock phase, is caused by the 
projectile compressing the fluid near the entry point to 
extremely high pressures. These pressures can range in the 
hundreds of thousands of pounds per square inch and cause a 
shock wave to move radially outward from the point of entry. 
However, the shock attenuates rapidly and, in general, is 
too weak to cause damage to the side and rear walls of the 
fuel tank; but it can cause outward bending of the wall in 
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the vicinity of the entry point and start cracks which can 
lead to eventual destruction of this wall. 

At any given point in the tank a (relatively) slow 
pressure pulse begins when the attenuated shock wave reaches 
that point and continues with the passage of the projectile 
through the liquid. This phase is called the pressure field 
phase. Its pressures range in the hundreds of pounds per 
square inch and last for several milliseconds. 

The third phase of hydraulic ram is called the cavity 
phase. As the. projectile moves through the fluid it must 
displace the fluid from its path. A large radial velocity 
is thus imparted to the fluid, and, as a result, the fluid 
moves away from the projectile path creating a cavity. 
Existing fluid pressures (gravity and pressures caused by 
tank boundaries) eventually halt cavity growth and cause the 
onset of collapse. Inside the cavity are fluid vapors and 
trapped air which entered during wall penetration. As the 
cavity collapses this air is compressed to a high pressure 
resulting in re-expansion of the cavity. Expansion and 
collapse repeat several times before all the stored energy 
is dissipated. Pressures during this phase are of the order 
of one thousand pounds per square inch. 

These three phases combine to impart a large velocity 
to the upper free surface of the liquid, if ullage is present, 
which in turn imparts an impulsive load on the top of the 
tank. This fourth phase of hydraulic ram is known as the 
free-surface phase . 
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The structural response to hydraulic ram Is brought 
about by the transfer of energy (projectile to fluid, fluid 
to walls, etc.). The failure mode is a function of how this 
energy is distributed among the four phases. In this study 
the projectile energy loss during entry wall penetration was 
investigated in detail and compared with experimental 
results. There are two areas of bullet penetration studies: 
those assuming a plate of finite thickness which the projec- 
tile passes through and those concerning a plate of. thick- 
ness greater than the penetration depth. The effect of 
projectile shape on penetration distance and energy loss by 
a projectile passing through a finite thickness plate were 
considered in this study. In order to be able to predict 
energy transferred to the fluid by the projectile during 
the shock phase and later during the cavity phase, a know- 
ledge of energy losses during entry wall penetration is 
essential. 

Knowledge of pressures during the three final phases of 
hydraulic ram is important because the major damage occurs 
during them. Fiber optic pressure transducers were designed 
in order to experimentally measure these pressures. The 
magnitude and duration of pressures during the shock and 
cavity phases were examined, and transducers were designed 
which are capable of accurately measuring these characteristics. 
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II. ANALYTICAL PENETRATION STUDIES 



In order to predict the phenomena taking place as a 
projectile impacts and passes through a plate, an analytical 
study of the entry wall phase of hydraulic ram is presented. 

An essential assumption of this study of projectile 

impact, and the assumption made by Dunn in Ref. 2, is that 

"there exists in each unconstrained material body at constant, 

uniform initial temperature a dynamic yield strength, a , 

* <y 

that is a constant of the material." Therefore there is one 
well-defined material constant which can be used to predict 
penetration data. 

This dynamic yield strength, when multiplied by the area 
of the projectile, A(x), which resists the motion, gives an 
estimate of the force retarding the bullet. An additional 
assumption is that the dynamic yield strength of the 
projectile is much greater than that of the target so that 
no substantial deformation of the projectile occurs. For 
these assumptions the retarding force is given by 



F = -o A(x) (1) 

y iji 

where a is the yield, strength of the target. Thus the 

y ip 

projectile equation of motion is 



x + 




A(x) = 0 



( 2 ) 



where x is the acceleration of the projectile mass m. One 
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salient feature of equation (2) is that the force retarding 
the projectile is only a function of x. It is assumed that 
there is no plastic deformation of the target material and 
no significant frictional force during penetration. 

Equation (2) was used for two types of penetration solu- 
tions: those concerning impact into a plate of thickness 

greater than penetration depth, and those concerning the 
passage of projectiles through thin plates- 

For the first case, multiplying (2) by x and integrating 
yielded 

P 

(^) 2 (* 2 - V, 2 ) = - / A(x) dx (3) 

u 1 0 



where G is the shear modulus of the plate and approximately 
equal to a /15 [Ref. 2], V is the impact velocity, and P 

y ip i 

the penetration depth. 

Equation (3) was solved numerically by a computer 
program. The program, for a given impact velocity, deter- 
mined the bullet trajectory parameter V as a function of 
incremental values of x until zero velocity was reached at 
the penetration depth. Figure 1 shows penetration depth vs. 
impact velocity for forty-five grain projectiles of three 
shapes. As expected, projectiles with a smaller cross- 
sectional area gradient penetrated further for a given impact 
velocity . 

The second penetration problem considered was that for a 
given plate thickness. Unlike the first case, where the 
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effective retarding force area distribution A(x) increased 
until it reached a constant maximum value, A(x) for a thin 
plate of given thickness increased only until the nose of the 
projectile penetrated the backside of the plate. At that 
time A(x) decreased due to projectile protrusion from the 
backside of the plate. Figure 2 shows two stages of pro- 
jectile penetration. It should be noted that the x- 
coordinate moves with the projectile so that x is defined 
as the distance from the nose of the bullet to the front 
side of the plate. In Figure 2a the bullet has partially 
penetrated the plate. The effective retarding area is shown 
as In Figure 2b the nose of the projectile has 

penetrated the back side of the plate and A g ^ is the 
annular area shown. Thus, as the maximum cross-sectional 
area reaches the backside of the plate, the area resisting 
motion goes to zero. It is at this point that penetration 
is complete. The minimum impact velocity which a projectile 
must possess to ensure its reaching this point is defined as 
the ballistic limit. 

The work done on the projectile during its deceleration 
in the plate is equal to the change in kinetic energy during 
penetration. An integration of the total retarding force 
during penetration yeilds 



W 




A 



max 



t 



(U) 



where t is the thickness of the plate. It is seen that the 
volume swept out by the projectile determines the penetration 
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energy loss, and the energy loss is totally independent of 

projectile shape. Alternately stated, this result predicts 

that, for a projectile of given mass and A , the ballistic 

max' 

limit is independent of projectile shape. 

With this in mind, equation (2) was solved assuming a 
constant cross-sectional area, that is, assuming a cylindri- 
cally shaped projectile. The equation of motion was inte- 
grated twice with appropriate initial conditions substituted 
to yield: 




T - 



GA t£ 
15m 2 



(5) 



where T is time and A is the projectile cross-sectional area. 

For a projectile with a velocity of impacting a plate 
of thickness t and exiting with a velocity of V , from (5) 
above , 



t 



V ± T 



GA 

15m 2 



( 6 ) 



and 




GA 

15m 



(7) 



Solving (7) for T and substituting into (6) and solving for 
V e yields 

v e ■ C v ! 2 - T?5 C ]1/2 (8) 



At the ballistic limit, the exit velocity is zero, hence 
equation (8) was solved for (denoted as V o to designate 
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ballistic limit): 



_ r GA . -.1/2 

o L 15m t J 



( 9 ) 



A relationship between the impact velocity, the exit 
velocity, and the ballistic limit can be found by substituting 
(9) into (8) and solving for V e : 




V i 2 

c (^> 2 

o 




( 10 ) 



Equation (10) shows, for impact velocities much higher 
than the ballistic limit, that the exit velocity approaches 
the impact velocity. 

Equation (10) may also be used to find an estimate of 
the percentage of energy loss for a projectile which pene- 
trates a finite thickness plate; 



* E lost “ [ 



E. 



] 100 = [ 



GA 

15m 



t ] [ 



] 100 ( 11 ) 



where E q is the projectile kinetic energy when traveling at 
the ballistic limit, and is the energy at impact velocity. 
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FIGURE 1. Penetration Depth vs. Impact Velocity for a Semi-Infinite Aluminum Target 
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FIGURE 2. Stages of Projectile Penetration 
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III. EXPERIMENTAL STUDY OF ENTRY WALL PENETRATION 



Figure 3 shows the basic components of the ballistic 
range used in the experiments. They were the rifle and its 
mount, chronograph screens and counters, aluminum test plate 
and holder, shadowgraph station, and bullet catcher. 

The rifle mounting system, shown in Figure 4, consisted 
of the rifle mount and an adjustable rigid rifle mount 
stand. The system was designed to hold rifles of .22, .222, 
and .30 caliber. In this study, a .222 Remington was used. 
The bullets utilized were of a mass of forty-five grains and 
of two shapes, hornet and spitzer. Projectile impact 
velocities ranged from 2^00 to 3000 feet per second. 

The velocity was measured by three chronograph screens, 
AVTRON No. A91^T333, mounted six feet apart. Each screen 
had a five-volt D.C. signal across it. When the screen was 
broken by a bullet, a pulse was sent to a counter. The first 
sensor started the first counter. The second sensor simul- 
taneously stopped the first counter and started the second. 
The third sensor stopped the second counter. Thus the 
average velocity between counters was obtained, allowing a 
prediction of the impact velocity at the plate. A typical 
chronograph station is shown in Figure 5- Also shown in 
this figure is a shadowgraph apparatus which was not utilized 
during this phase of testing. The counters used were 
Monsanto 101B 1MHz. The trigger error was less than i 0.3? 
of one period for sine waves with signal-to-noise ratio of 
40 db or better. 
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The empty fuel tank entry wall was simulated by a three- 
inch by four- inch 7075-T6 aluminum plate. Three thicknesses, 
0.05 inches, 0.09 inches, and 0.16 inches were tested. 

Early in the test phase it was determined experimentally 
that the size and method of supporting the test plate had no 
significant effect on the exit velocity of the projectile. 

In addition to stopping the second counter, sensor three 
was used to trigger the delay unit which generated a pulse 
which could be time delayed from 190 yseconds to 1700 
yseconds. The delay unit circuitry is shown in Figure 6. 

The delayed pulse was then used to trigger a spark source 
just after bullet penetration of the plate.. The light from 
the spark source was directed to a film plate by a concave 
mirror on the opposite wall. The mirror collimated the light 
which then passed behind the test plate and onto the film. 

In this manner, shadowgraphs were taken to examine the shape 
of the bullet after it had passed through the plate. The 
shadowgraphs also showed the resulting flow field behind the 
plate. Figure 7 shows the spark source and mirror as they 
were set up. 

The chronograph screens mounted directly behind the test 
plate were used to measure the velocity of the bullet as it 
exited the plate. Figure 8 is a photograph of this apparatus 

The bullet catcher is shown in Figure 9- The sides, top, 
and bottom were 1/2-inch aluminum. The three-foot square 
front had three 3/8-inch plywood baffles Inserted into it. 

The back of the catcher was made of two plates of 3/8-inch 
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steel mounted at an angle of 45 degrees to the path of the 
bullet. This ensured bullet deflection downward into a 
layer of sand, reducing the likelihood of ricochets. 
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FIGURE 3. Ballistic range components. 




FIGURE 4. Ballistic Range Rifle Mounting System 



20 





#■ 




FIGURE 5- Ballistic Range Chronograph Station 
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FIGURE 6. Delayed Pulse Generator Circuitry (Variable) 
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FIGURE 7- Spark Source and Mirror 
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FIGURE 8. Test Plate and Chronograph Screens 
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FIGURE 9- Ballistic Range Bullet Catcher 



IV. FIBER OPTIC PRESSURE TRANSDUCER DESIGN 



A. INTRODUCTION 

An accurate means of measuring pressures generated 
during all phases of hydraulic ram is needed to assess the 
validity of any analytical studies. The most popular types 
of pressure measuring devices in use today are the semi- 
conductor strain gauge and the piezoelectric transducer. 

Both costly and fragile, these devices are also sensitive 
to temperature changes and contamination from small dust 
particles. A fiber optic pressure transducer was designed 
in order to overcome these problems. Two optical glass 
fiber bundles were used, one to transmit light to the inner 
surface of the transducer diaphragm and one to collect and 
transmit reflected light to a photomultiplier tube. This 
signal output was then amplified and could be measured from 
a photograph of the resulting oscilloscope trace . The 
intensity of the output light, being a function of displace- 
ment between the fiber bundle and the diaphragm, was thus 
made a measure of pressure on the diaphragm. A possible 
arrangement of the pressure measuring apparatus is shown in 
Figure 10. 

The fiber optic transducer is in reality a proximity 
gauge. Because pressure on the diaphragm causes it to 
deform slightly and decrease the distance between the 
diaphragm surface and the fiber ends, output can be read in 
terms of pressure increase rather than distance decrease. 
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Figure 11a, obtained from Ref. 6, shows two fibers, one 
emitting light and one collecting light. As the displacement 
between the fiber ends and the diaphragm increases, the cone 
of light from the transmitting element illuminates an 
increasingly larger area of the diaphragm. The light 
reflected from the diaphragm becomes a cone of light, 
increasingly intense as the diaphragm moves away from the 
fiber ends . Figure lib shows typical outpuit of the 
receiving fibers vs. fiber-diaphragm displacement. The 
front side of the curve indicates linearity between output 
and distance. As the distance increases, however, the cone 
of light from the transmitting fiber fills fche entire 
diaphragm, and the curve reaches its peak. Beyond this 
point, the intensity of light to the receiving element 
decreases as the inverse of the distance squared [Ref. 6]. 

B. DESIGN CONSIDERATIONS 

Due to the great disparities between the peak pressures 
and the time scales of the pressure signals :In the shock and 
cavity phases of hydraulic ram, two separate pressure trans- 
ducer designs are necessary. Figure 12 shows the basic 
geometry of both transducers. In addition ho the optical 
fiber probe, the transducers consist of an outside tube and 
a diaphragm silver soldered or welded to the tube end. A 
sleeve is required inside the tubing to ensure orthogonal 
orientation of the probe tip to the diaphragm. 
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Linear response of the pressure transducers is desired, 
and in order to achieve this the maximum deflection of the' 
diaphragm is limited to less than half its thickness. This 
allows the light transmission to remain on the linear 
portion of the displacement curve of Figure 11. The resonant 
frequency of the diaphragm must be large enough to provide 
adequate frequency response to the transducer. Finally, the 
pressure transducer must be able to withstand the maximum 
pressures occurring during hydraulic ram and should be placed 
so that the diaphragm is flush with the tank wall when 
measuring shock phase pressures. The transducer should be 
isolated from wall vibration, if any; therefore, it was 
insulated by a seal made of RTV. 

C. SHOCK PHASE DESIGN 

Figure 13 shows peak pressures (predicted by the strong 
shock theory of [4]) encountered at various distances from 
point of impact of a .30 caliber projectile penetrating the 
fluid at 3000 fps and an energy of 34,000 ft-lb. Measurable 
pressures are of duration of about 12 ysec and a pressure 
transducer placed approximately 1.4 inches from the point of 
initial impact will experience a peak pressure of 7.8 x 10 1 * 
psi, 12 ysec after impact. Figure 13a indicates that the 
shock pressure signature can be approximated by the sawtooth 
function shown in Figure 13b. The Fourier Expansion of the 
sawtooth is given by 



28 



f(x) = 



1 

TT 



n 



k=l,3,5..Tr 2 k 2 



cos kx + 



(-ly 

irk 



sin kx 



} 0 <x< 



0<x<2tt (12) 



where a k 



"272 311(3 

ir k 



(-D k 

TTk ' 



For adequate frequency response of the transducer, it 
is desired that the amplitude of the high frequency signal 

content be at most ten percent of the fundamental amplitudes 

_2 

a k and b k * For a k proportional to k , the smallest value of 
k that meets this requirement is 5* However, b k is propor- 
tional to k"* 3- and thus the smallest value of k is 11. The 
high frequency content of the shock pressure signal is then 
approximately k/24 ysec or 0.46 MHz. 

Requiring the transducer to be capable of measuring peak 
pressure of 100,000 psi. Figure 15 , obtained from Ref. 3, 
shows that the minimum diaphragm thickness is 23 mils. 

According to Figure 16 , all stainless steel diaphragms of 
sufficient diameter (greater than 0.160 inch to accommodate 
the . 109 -inch fiber bundle) have resonant frequencies near 
the operating frequency. A diaphragm 35 mils thick and 0.160 
inches in diameter has a resonant frequency of 0.5 MHz and 
has a maximum deflection of just under 1 mil at 100,000 psi. 
Because the fiber optic system output resolution is approxi- 
mately 0.01 mils, only pressures above 1000 psi can be measured 
for this proposed design. 
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D. CAVITY PHASE DESIGN 



The geometry of a transducer measuring pressures during 
the cavity phase must be different from that required to 
measure pressures during the shock phase due to lower pres- 
sures acting over longer periods of time during cavity 
collapse . 

Figure l4a, obtained from Ref. 5> indicates, for a 0.30 
caliber projectile, maximum cavity collapse pressures were 
less than 1000 psi and each pulse lasted for approximately 
0.5 msec. The pressures rose and fell at nearly the same 
rate, allowing the pressure signature to be approximated by 
a triangular wave, as shown in Figure i4b. 

As in the previous analysis, the Fourier expansion of 
the triangular function yields: 

“J h ft -1 

f(x) - -k - — 5 - \ — p cos kx 0 <x<_2tt (13) 

* TT k=l, 3, 5 ... k 

Again the amplitude of the high frequency components must be 
10 # of the fundamental amplitude for adequate frequency 
response. The frequency response required by this pressure 
signal is then k/0.5 msec or .01 MHz. 

A stainless steel diaphragm 10 mils thick with a 
diameter of 0.190 inches has a resonant frequency of very 
near 0.1 MHz, and is capable of withstanding peak pressures 
of 3000 psi. It is suggested that this geometry be used 
for the cavity phase transducer. 
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FIGURE 10. Experimental Arrangement of Fiber Optic Pressure Transducer 
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FIGURE 11. Displacement Measurement 
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FIGURE 12. Transducer Cross Section 
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FIGURE 13. Pressure Pulses and Sawtooth Representation 
for Shock Phase Pressure Signal 
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FIGURE lU . Diaphragm Maximum Deflection 
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FIGURE 15. Diaphragm Maximum Deflection 
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FIGURE 16. 



Diaphragm Fundamental Resonant Frequency 



V. RESULTS AND CONCLUSIONS 



One important assumption made in the analytical study was 
that the yield strength of the projectile was significantly 
higher than that of the target material. However, the pro- 
jectile was copper jacketed, and the yield strength of copper 
is approximately 0.39 x 10^ psi as compared with 0.25 x 10^ 
psi for aluminum. Figure 17 shows clearly that, although 
the projectile and target yield strengths are of the same 
order of magnitude, projectile shape was not altered during 
plate penetration, except very close to the nose of the 
projectile. Thus it can be concluded that for hypervelocity 
impact phenomena, the assumption only need be made that the 
yield strengths of the projectile and target materials be 
of the same order of magnitude. 

The analytical result that projectile shape is not 
material in predicting penetration energy losses is also 
shown in the experimental results. Equation (9) was derived 
using that assumption. Figure 18 shows a plot of exit veloc- 
ity vs. impact velocity for 45-grain projectiles of two 
shapes penetrating aluminum plates of three thicknesses. 

The solid lines are plots of Equation (10), using Equation 
(9) as the definition for ballistic limit. Because both 
shapes yielded similar data, no distinction between them was 
made. Data points agreed well with prediction leading to 
the conclusion that projectile shape does not, in fact, 
significantly affect energy loss through a plate. 
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The theoretical energy loss from Equation (11) was 
plotted vs. impact velocity in Figure 19. Experimental 

I 

data again agree closely with prediction, reinforcing the 
concept that projectile shape has no significant effects 
on energy loss through thin plates. 

It can be concluded, therefore, that a means is avail- 
able whereby energy losses through a wall of an unfilled 
fuel tank can be predicted accurately and simply. This is 
the first step in finding energy transferred during the 
hydraulic ram phenomenon. Effects of fluid foundation on 
the penetration characteristics of the walls are expected 
to be minimal, but further experimental study of this effect 
should be undertaken. 

It should be noted that all experimentation and analyti- 
cal work in this study has been conducted using bodies of 
revolution that do not tumble. However, an application of 
the concept of the shape independence of energy loss can be 
made in the case of fragment penetration. This representa- 
tion will only be applicable to early stages of hydraulic 
ram, as projectile shape does have major effects on the size 
and intensity of the cavity phase . Aircraft damage and 
penetration characteristics of fragments (such as from SAM 
warheads) may be approximated by cylindrical bodies of equal 
mass and maximum cross-sectional area. It is expected that 
fragment energy loss during wall penetration will follow 
Equation (11). 
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FIGURE 17. 



^5 Grain Spitzer Entering and Exiting 
Aluminum Plate 
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VI. RECOMMENDATIONS FOR FURTHER STUDY 



Experimental studies in the area of hydraulic ram are 
continuing and further information will be gained. A few 
specific recommendations are presented as guidelines for 
f ut ure research. 

The ballistic range is equipped to fire .30 caliber 
projectiles. This should be done in order to obtain data 
for projectiles of higher energies and larger damage radii. 
In addition, .30 caliber projectiles offer a greater 
variety of nose shapes which can add information to the 
previous data on bullet shape effects. 

Filled tanks should be penetrated by projectiles of 
various caliber and mass. This will yield information 
concerning the effect of fluid foundation on penetration 
of thin walls: 

Finally, it is recommended that the pressure transducers 
be installed and that experiments be conducted to confirm 
their effectiveness. If successful, this would give 
pressures as functions of time for the shock and cavity 
phases of hydraulic ram, a major step in understanding the 
phenomenon. 
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APPENDIX A 



Explanation of Computer Program 

This program computes the distance which a projectile of 
arbitrary shape and mass will penetrate into a semi-infinite 
metal plate. Inputs to this program are the shear modulus 
of the plate material, the mass of the bullet, its maximum 
radius, and its initial velocity. The cross-sectional area 
of the bullet at each .032 inch aft of the nose is also read 
in. Equation (3) is integrated trapezoidally using the 
bullet area distribution and velocity is tested until it 
becomes negative. The previous velocity is essentially zero, 
and the corresponding distance is output as the penetration 
depth. 
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C ****************************************** *************** 
C * 

C * THIS PROGRAM WILL YIELD FOR A GIVEN INITIAL VELOCITY 

C * THE DISTANCE THAT A PROJECTILE OF ARBITRARY AREA 

C * DISTRIBUTION WILL PENETRATE INTO A PLATE CF 7075-T6 

C * ALUMINUM 

C * 

c **************************************** ******q(* ********* 
DIMENSION AREA(300G) , A INTEG ( 3000 ) ,X(3000) ,A(15) , 

ID I ST( 15 ) 

READ IN INITIAL DATA 

DIST( 1) =0. 

RAC=. 112 
L=2999 
ALENG=. 416 
N = 14 

AINTEGI 1) =0. 

DELX= 1. E-4 
C=15 . 1 1 E09 
LPL=L+ 1 
DO 40 I =2 * N 
40 DIST(I) = (I-1)*.032 

READ IN AREA DISTRIBUTION 

READ<5, 10) ( A ( J) , J=1,N) 

10 FORMAT ( SF10. 8 ) 

DO 100 1=1, LPL 
X(I)=DcLXMI-l) 

I F < X < I) .GT.ALENG) GO TO 99 
GO TO 88 

99 AREA! I ) =3. 141 59* ( RAD** 2 ) 

GO TO 100 

PIF2 IS AN INTERPOLATION FORMULA 

88 AREA! I) =P I F2 ( X, DI ST , N , A , I ) 

100 CONTINUE 
K = I 

READ IN GIVEN IMPACT VELOCITY 

4CC READ! 5 » 800) VINIT 
800 FORMAT (F 8.1) 

IF(VINIT. GT. 36000. ) GO TO 850 
WRITE (6*3 00) 

300 FORMAT < 'ODISTANCE INTO PLATE ' , 12X , • VELOC I T Y • ) 

WRITE (6 *350) VINIT 
3 50 FORMAT ( • 0 • , 1 1 X , ' 0. 0 • , 13X , F 1 5. 4 ) 

NUMERICAL INTEGRATION 

DO 600 1=1, K 

AINTEGt 1+1 )=AIMTEG(I )+( ( AREA( I + 1 ) +AREA( I ) )/2. )*DELX 
E = V I N I T * ? 2~ ( 2 . * C * A I NT E G ( 1*1) ) 

IF(E.LT .0. ) GO TO 200 
6 OC CONTINUE 
200 E = VINIT-1'2~(2.*C*AINTEG( I ) ) 

V EL=S GRT ( E) 

WRI TE ( 6 ,550) X(I),VEL 
550 FORMAT ( '0'*2X,F15.4,10X,F15.4) 

GC TO 400 
85C STOP 
END 
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FUNCTION PIF2(X,XLIST,N,FLIST,J) 

SUBPROGRAM FUNCTION 2ND ORDER LOOKUP 
DIMENSION XLI ST( 200) ,F LI ST (200) ,X(2Q0) 
BLIF(P»Q,R»S»T)=( (Q-P)*(S-T)/(R-Q)+S) 

I F ( X( J ) —X LIST (N) ) 2,1,1 

1 I=N— 1 
GC TO 5 

2 I F ( X( J) —XLI ST C 1 ) ) 4,4,6 

10 IF (K— 1 ) 11,11,12 

4 1=1 

5 K = 1 

GO TO 30 

6 K=2 

7 DO 8 1=1, N 

I F ( X( J ) -XL 1ST ( I ) ) 9,y ,8 

8 CONTINUE 
I =N 

9 1=1-1 

3C BLIF1=BLIF(X( J) , X LI ST ( II , X L 1ST 1 1 +U , FL I ST ( I), 

1FL 1ST ( 1 + 1) ) 

11 PIF2 = 3LIF1 
RETURN 

12 I F ( ( 1 + 2 )— N ) 13,13,16 

13 I F ( ( I— 1 ) — 1) 15,14,14 

14 IF(ABS(XLIST(I-1)-X(J) ) -ABS (XL I ST ( 1+2 )-X ( J ) )) 16,15,15 

15 L = 1+2 

GC TO 17 

16 L=I-1 

17 BLIF2=BLIF(X( J) , XL IST(I) ,XLIST( L) ,FLIST(I) ,FLIST(L)) 
PIF2=BLIF(X(J),XLIST( 1 + 1) ,XLIST( L ) ,BL IF 1 , BL IF 2) 

18 RETURN 
END 
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